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R.H.W. HOPPE, W.G. LITVINOV. 



Abstract. We set up and study a coupled problem on stationary non-isothermal flow of 
electrorheological fluids. The problem consist in finding functions of velocity, pressure and 
temperature which satisfy the motion equations, the condition of incompressibility, the 
equation of the balance of thermal energy and boundary conditions. 

We introduce the notions of a P-generalized solution and generalized solution of the 
coupled problem. In case of the P-generalized solution the dissipation of energy is defined 
by the regularized velocity field, which leads to a nonlocal model. 

Under weak conditions, we prove the existence of the P -generalized solution of the 
coupled problem. The existence of the generalized solution is proved under the conditions 
on smoothness of the boundary and on smallness of the data of the problem 



1. Introduction 

Electrorheological fluids are smart materials which are concentrated suspensions of polar- 
izable particles in a nonconducting dielectric liquid. In moderately large electric fields, the 
particles form chains along the field lines, and these chains then aggregate to form columns 
[11]. These chainlike and columnar structures cause dramatic changes in the rheological 
properties of the suspensions. The fluids become anisotropic, the apparent viscosity (the 
resistance to flow) in the direction orthogonal to the direction of electric field abruptly in- 
creases, while the apparent viscosity in the direction of the electric field changes not so 
drastically. 

The chainlike and columnar structures are destroyed under the action of large stresses, 
and then the apparent viscosity of the fluid decreases and the fluid becomes less anisotropic. 

On the basis of experimental results the following constitutive equation was developed in 
[5]: 

aij{p,u,E) = -p6ij + 2ip{I{u), \E\,ij,{u,E))eij{u), i,j = l,...,n, n = 2 or 3. 

(1.1) 

Here, aij{p,u, E) are the components of the stress tensor which depend on the pressure p, 
the velocity vector u = (ui, . . . , u„) and the electric field strength E = [Ei, . . . , En), Sij are 
the components of the unit tensor (the Kronecker delta), and Sij{u) are the components of 
the rate of strain tensor 

1 , dui duj s 
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Moreover, I{u) is the second invariant of the rate of strain tensor 

n 

I{u)=J2i^iji^)f, (1-3) 
and ip the viscosity function depending on I{u), \E\ and fj,{u,E), where 

So fi{u, E) is the square of the scalar product of the unit vectors |^ and -j^. The function 

is defined by (1.4) in the case of an immovable frame of reference. If the frame of reference 
moves uniformly with a constant velocity u = (ui, . . . , -tin), then we set: 

As the scalar product of two vectors is independent of the frame of reference, the constitutive 
equation (1.1) is invariant with respect to the group of Galilei transformations of the frame 
of reference that are represented as a product of time-independent translations, rotations 
and uniform motions. 

It is obvious that fi{u, E){x) G [0, 1], and for fixed yi, y2 £ where M+ = {z G M, z> 
0}, the function j/3 (^(yi, 1/2, ys) reaches its maximum at ys = and its minimum at 
y3 = 1 when the vectors u{x) + u and E arc correspondingly orthogonal and parallel. 

The function fi defined by (1.4), (1.5) is not specified aX E = Q and at tt = 0, and there 
does not exist an extension of /x by continuity to the values of u = and E = {). However, 
at = there is no influence of the electric field. Therefore, 

</f(yi,0,y3) = <p(yi), y3e[0,l], (1.6) 

and the function /i(u, E) need not be specified at = 0. Likewise, in case that the measure 
of the set of points x at which u{x) = is zero, the function jjl need not also be specified at 
n = 0. But in the general the function is defined as follows: 

( ai + u{x)+u E{x) 

where / denotes a vector with components equal to one, and a is a small positive constant. 
If u{x) ^ almost everywhere in Q,, we may choose a = 0. 

The viscosity function is identified by approximation of flow curves, see [5], and it was 
shown in [5] that it can be represented as follows: 

^{I{u), \Elii{u, E)) = h{\Elii{u, E)){\ + I{u))-h + i,{I{u), \E\,ti{u, E)), 

(1.8) 

where A is a small parameter, A > 0. 

In the special case that the direction of the velocity vector u{x) at each point x at which 
E{x) 7^ is known, the function x — iJ,{u,E){x) becomes well-known, and the viscosity 
functions (1.8) takes the form 

ip{I{u), \E\,x) = e{\E\,x){X + 7(«))-5 + Vi(/(n), \E\,x), (1.9) 
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where 



e{\E\,x)=b{\E\,f,{u,E){x)), 

V'l (/(«), 1^1, x) = ij{I{u), \E\,fi{u, E){x)). 



(1.10) 



Here and thereafter the function of electric field E is assumed to be known. The equations 
for the functions E and {p, v) arc separated, see [5]. 

In actual practice, the temperature fields in electrorheological fluids are nonuniform, and 
in many cases this non-homogeneity drastically affects on the flow of electrorheological fluids. 
Because of this, we reckon below that the functions b, if), e, ■^i in (1.8) and (1.9) depend also 
on the temperature r. Therefore, we consider the following viscosity functions 



Further we study coupled problems on non isothermal flow of electrorheological fluids with 
the constitutive equation (1.13) for k equals one and two. 

In Sections 2 and 3 we introduce equations governing the process of non isothermal flows 
of electrorheological fluids. The coupled problem on non-isothermal flow consists in finding 
functions of velocity, pressure and temperature which satisfy the motion equations, the 
condition of incompressibility, the equation of the balance of thermal energy and boundary 
conditions. 

We introduce the notions of a P-generalized and generalized solution of the coupled 
problem. In case of the P-generalized solution the dissipation of energy is defined by the 
regularized velocity field, which leads to a nonlocal model. 

Some auxiliary results are set in Section 4, and the existence of a P-generalized solution 
of the coupled problem for the viscosity function is proved in Section 5. 

Assuming that the data of the problem are small and the boundary of the domain under 
consideration is smooth, we prove in Section 6 the existence of a generalized solution of the 
coupled problem for the viscosity function (^2- 

In Section 7 we consider coupled problem for the viscosity function ipi and formulate a 
result on the existence of the P-generalized solution for ip\ . A proof of this result is contained 
in Sections 8 and 9. 

2. Equations of flow and thermal balance, boundary conditions. 

We consider stationary problem on non-isothermal flow of the electrorheological fluid. 
The inertia forces are assumed to be small as compared with the internal forces caused by 
the viscous stresses. Then the equations of motion take the following form: 

- 2—\ipkeij{u)\ = K, in 17, i = 1, . . . , n, fc = 1 or 2, (2.1) 

where Ki are the components of the volume force vector K. In (2.1) and below the Einstein 
convention on summation over repeated index is applied. 

The velocity function u meets the incompressibility condition 



^1 = h{\E\,^l{u, E), r)(A + /(«))-^ + ^{I{u), \E\,fi{u, E), r), 
V2 = ei\E\, T, x) (A + I{u)y'^ + (/(n) ,\E\,t,x), 



(1.11) 
(1.12) 



and the stress tensor is defined by 

aij = -pSij + 2(pk£ij{u), A; = 1 or 2, i,j = l,...,n. 



(1.13) 




n 



(2.2) 
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The equation of the balance of thermal energy is the following: 

dx'i 



X2_/^ +^^'^'^^('") " ^ in J^, A; = 1 or 2. 

1=1 " ' 



(2.3) 



Here X is the thermal diffusivity, e = {pcp)~^, where p is the density, Cp specific heat. We 
reckon that x ^-nd e are positive constants. 

The first, second and third terms in the left-hand side of (2.3) determine the increments 
of the temperature in the unit of time produced respectively by heat conduction, dissipation 
of energy and thermal convection. 

We assume that is a bounded domain in M", n = 2 or 3 with a Lipschitz continuous 
boundary S. Suppose that and are open subsets of S such that Si is non-empty, 
S*! n 5*2 = 0, and 5*1 U 5*2 = S. We consider mixed boundary conditions for the functions 
u,p, wherein velocities are specified on and surface forces are given on ^2, i.e. 



u 



Si 



u, 



-pSij + ipkeij{u)]i'j 



Si 



1, . . . , n, = 1 or 2. 



(2.4) 
(2.5) 



Here Fi and Uj are components of the vectors of surface force F = (i^i, . . . , i^„) and the unit 
outward normal = (i^i, . . . , f„) to 5, respectively. 

Temperature field on the boundary S is considered to be given 



T. 



(2.6) 



We assume that u E H-2 (Si), t E H-2 (S). Then there exist functions u and f such that 



u 



u £ i?^(0)". 



Si 



divu = 0, 



= r. 



Suppose also 



(2.7) 
(2.8) 

(2.9) 



3. P-GENERALIZED SOLUTION OF THE COUPLED PROBLEM FOR THE FUNCTION (f2, AND 

THE EXISTENCE THEOREM. 

We consider coupled problem for the function (p2 defined by (1.12). We assume that the 
functions e and ipi satisfy the following conditions: 

(CI): e : {yi,y2,x) — e{yi,y2,x) is a function continuous in M^. x R x $7 and in 
addition 

0<e{yi,y2,x) <ao, {yi,y2,x) e R+ x R xTl, (3.1) 

oo being a positive constant. 

(C2): "01 : yi,y2, 2/3, x "01 (yi , 2/2 , 2/3 j x) is a function continuous in ]R_|_ x M_|_ x R x J7, 
and for an arbitrary fixed {y2, y3,x) G R+ x R x $7 the function il;i{.,y2,y3, x) : yi 
fpi{yi,y2,y3, x) is continuously differentiable in R_|_, and the following inequalities 
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hold: 

02 > V'i(yi,2/2,y3,a;) > ai, (3.2) 

dip 

iki{yi,y2,y3,x) + 2— {yi,y2,y3,x)yi > 03, (3.3) 
■^(2/1,2/2,2/3,3;) yi < 04, (3.4) 

where Oj, 1 < z < 4, are positive constants. 

The inequalities (3.1) and (3.2) indicate that the viscosity function is bounded from below 
and above by positive constants. (3.3) implies that for fixed values oi\E\, t, x the derivative 
of the function I{u) — > G{u) is positive, where G(u) is the second invariant of the stress 
deviator G{u) = 4ip2 I{u). This means that in the case of simple shear flow the shear stress 
increases with increasing shear rate. (3.4) is a restriction on for large values of yi. The 
inequalities (3.2)~(3.4) are natural from the physical point of view. 

We consider the following spaces: 

X = {u\u e H^in)"", u\s, = 0}, (3.5) 
V = {u\ueX, div'u = 0}. (3.6) 

By means of Korn's inequality, the expression 

\\u\\x= ( / I{u)dxy (3.7) 

defines a norm on X and V being equivalent to the norm of H^[^)"'. 

Everywhere below we use the following notations: if y is a normed space, we denote by 
Y* the dual of Y, and by (/, h) the duality between Y* and Y, where / G Y*, h gY. In 
particular, if / G L2{^) or f E L2(il)", then (f-h) is the scalar product in L2{^) or in 
L2{^)^, respectively. The sign^ denotes weak convergence in a Banach space. 

Define a bilinear form tt in H^{Q) x H^{0,) as follows: 



dOi 002 

•^^ i=l 

The expression 



<Gi^02) = J^Y.'^^dx, ei,e2eH\n). (3.8) 



\\9\\i = {7T{e,e))-^ (3.9) 

defines a norm in Hq^Q) that is equivalent to the norm of H^(Q.). 
We introduce a mapping N : X x Hq (O) X* by 

{N{v,0,h) = 2 [ [e{\E\,f + C,x){X + I{u + v)y'^ 
Jn 

+'4>i{I{u + v),\E\,f + C,x)]eij{u + v)£ij{h)dx, v,heX, ( £ H^{n), 

(3.10) 



and we assume that A > in (3.10). 
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Determine an operator A : X x Hq{Q,) H ^(0), where H ^{Q) is the dual space of 
Hq{Q,), as foUows: 

(^A{v, 0, = ^(f + Oin + Vi)^ dx + 2ejje{\E\,f + C, x){X + I{u + v))-^ 

+Ml{i^ + v),\E\,f + C,x)]IiP{u + v))^dx-xAr,0, vex, {Q,i)eHl{nf. 

(3.11) 

Here P is an operator of regularization given by 

Pu{x)= Lo{\x - x'\)u{x')dx' , xen, (3.12) 

where 

a;GC°°(M+), supp w = [0, a], uj{z)>0 z eR+, 
/ a;(|x|)(ix = 1, a is a small positive constant. (3.13) 

In (3.12) we assume that the function u is extended to W^. 
We denote by B the operator div, i.e. 

Bu = dwu. (3.14) 

It is obvious that B G C{X,L2{0,)), and we denote by B* the adjoint to B operator. 
Consider the following problem: find a triple {v,p, (") such that 

{v,p,C) e X X L2{n) X H^{n), (3.15) 

(n{v, 0, h) - [b*p, h) = {K + F,h), hex, (3.16) 

(Bv,q)=Q, qeL2{n), (3.17) 

7r(C,0 - ^{Mv,0,C) =0, Ce H'oin). (3.18) 
Here we use the notations 

(K,h^ = J Kihidx, [F,h^ = J Fihids, heX. (3.19) 

The triple {u = u + v, p, r = f + C), where v,p, ^ is a solution of the problem (3.15)-(3.18), 
will be called the P-generalized solution of the problem (2.1)-(2.6) for the viscosity function 
(p2; in the special case that P is a unit operator the triple {u = u + v, p, t = f + C) is a 
generalized solution of the problem (2.1)-(2.6). 

Indeed, by use of Green's formula it can be seen that smooth generalized solution of the 
problem (2.1)-(2.6) is a classical solution of this problem. On the contrary, if (u,p,9) is a 
classical solution of the problem (2.1)-(2.6), then the triple {v = u — u, p, ( = r — f) is a 
solution of the problem (3.15)-(3.18), wherein P is the unit operator. 

From the physical point of view the presence of the operator P in (3.18) means that the 
value of dissipation of energy at a point x depends on the values of the rate of strain tensor 
at points belonging to some small vicinity of the point x. Therefore, the presence of the 
operator P in (3.18) is natural from the physical point of view and it leads to a nonlocal 
model. 
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Theorem 3.1. Let ^ be a bounded domain in M", n = 2 or 3 with a Lipschitz continuous 
boundary S, and suppose that the conditions (CI), (C2), (2.7) -(2.9) are satisfied. Then 
there exists a solution of the problem (3.15) -(3.18). 

4. Auxiliary results. 

The next lemma follows from the results of the work [5] . 

Lemma 4.1. Assume that the conditions (CI), (C2), (2.7), (2.8) are satisfied. Then for an 
arbitrary fixed ( G Hq(Q) the partial function N{.,Q : v — N{v,Q, where N{v,Q is given 
by (3.10), is a Lipschitz-continuous strictly monotone and coercive mapping of X into X* , 
i.e. 

\\N{v,0-N{w,0\\x* <^^l\\v-w\\x, v,weX, (4.1) 
{N{v,C) - N{w,C),v - w) > fi2\\v - w\\x, v,weX, (4.2) 
{N{v,C),v)> H3\\vfx - f^4v\\x, vex, (4.3) 

where 

fii = 2a2 + 4(04 + aoA~2), fi2 = min(2ai, 203), 

113 = 2ai, ii4 = 2{aoX~^ +a2)(^J I{u)dxy . (4.4) 

Lemma 4.2. Let Q be a bounded domain in W^, n = 2 or 3 with a Lipschitz continuous 
boundary S, and let the operator B G C{X, L2{^i)) be defined by (3.14) Then, the inf-sup 
condition 

inf sup „ i/^^^'if^ >Pi>0 (4.5) 

neL2{n) vex \\v\\x ||/^||L2(n) 

holds true. The operator B is an isomorphism from V'^ onto L2(f2), where V'^ is orthogonal 
complement ofV in X, and the operator B* that is adjoint to B, is an isomorphism from 
^2(0) onto the polar set 

1^0 = {/GX*, (/,«) = 0, ueV}. (4.6) 

Moreover, 

\\B-'\\ciL2in),v-)<^^, (4.7) 

ll(^*)"'ll£(vo,L.(n))<^- (4.8) 

For a proof see [1]. Lemma 4.2 is a generalization of the inf-sup condition in case that 
the operator div acts in the subspace ifQ(r2)" (see [3]). This result was first established in 
an equivalent form by Ladyzhenskaya and Solonnikov in [6]. 

Theorem 4.1. Suppose the conditions (CI), (C2), (2.7)-(2.9) are satisfied. Then for an 



arbitrary ^ € Hq{^) there exists a unique pair {v{Q, p{C)) such that 

viO G X, piO G L2{n), (4.9) 

{N{v{0,0,h) - {B*p{0,h) = {K + F,h), hex, (4.10) 

{Bv{O,q)=0, qeL2{n). (4.11) 



In addition 

\HC,)\\x< — {\\K + F\\v*+iia), C^Hl{n), (4.12) 

and the function C, — ^ (^^(C); P(C)) o, compact mapping of 11^(0,) into X x ^2(0). 
The condition 

a^Co inH^in) 

implies 

viCk) v{Co) in X, 

Proof. It follows from (4.10), (4.11) and (3.6) that the function v{Q is a solution of the 
following problem: 

viOeV, {N{v{0,0,h) = {K + F,h), heV. (4.13) 

Lemma 4.1 and the results on solvability of equations with monotone operators, see e.g. [7] 
imply that there exists a unique solution of the problem (4.13). Equality (4.13) yields 

N{v{C),0-K-FeV''. (4.14) 

By virtue of Lemma 4.2 and (4.14) there exists a unique function p(Q G L2{^) such that 
the pair (f (C), PiC)) is a solution of the problem (4.9)-(4.11). 
Taking h = v{C) in (4.10), we obtain 

{N{v{0,C),v{C)) = {K + F,v{C)) < \\K + F\\v*\\v{C)\\x, (4.15) 

and (4.12) follows from (4.3) and (4.15). Inequality (4.12) implies 

lb(C)l|L,(n)<c, CeH^oi^). (4.16) 

Let now 

Ck ^ Co in H^{n). (4.17) 

We take the following notations: 

v"' = v{Cm), Im = Iiu + vn, £^ = Si,' (ti + O , m = 0, 1, 2, . . . , 
ifmk = ei\E\,f + Ck, x)iX + I{u + ?;"*))-§ 

+V;i(/(tx + i;"^),|S|,f + a,rr), A;, m = 0, 1, 2, . . . . (4.18) 

It follows from (4.13) that 

{N{v^,(:k),h) = {K + F,h), heV, k = l,2,... 
{N{v''),Co),h) = {K + F,h), heV. 

Taking h = v'' — v^, we obtain from here by (3.10) that 

{^kk 4 - ^00 4)(4 - 4) dx = 0. (4.19) 
(4.19) implies 

+ = 0, (4.20) 



I' 



where 

Mf = / {^kk 4 - ^Ok 4)(4 - 4) dx, (4.21) 
Jn 

M| = / i^ok - Voo) 4(4 - 4) dx. (4.22) 
The Cauchy inequahty and the notation (4.18) imply 

|M|| < {J^i^ok - moflo d^yij^ ^(^' - dx) " . (4.23) 

By virtue of (4.17), we can extract a subsequence {Cq} from the sequence {(k} such that 

C ^ Co a.e. in 0,. (4.24) 
Now the continuity of the functions e and V'l (see (CI), (C2)) yields 

y^Or] Voo a.e. in Q. 

The functions iporj are bounded by a constant aoA~2 +02. Therefore, the Lebesque theorem 
and (4.23) imply that 

limM2' = 0. (4.25) 

By analogy with the stated above, we can see that from any subsequence {M|"} extracted 
from the sequence {M|}, one can extract a subsequence {M™'} such that limM^' = 0. 
Therefore, limM|^ = 0, and by (4.20), we get 

limMf = 0. (4.26) 

It follows from (3.10), (4.2), (4.18) and (4.21) that 



^ 2 



(iv(^;^a)-iv(^;^a),^ -^") 

This inequality together with (4.26) yields 

v{Ck)^v{Co) inX (4.27) 

We have 

\\N{v{Ck)Xk)-N{v{Co)Xo)\\x* 

< \\N{v{Ck),Ck) - N{v{Co),Ck)\\x* + \\N{v{Co)Xk) - N{v{Co),Co)\\x*. 

By virtue of (4.1) and (4.27) the first term on the right-hand side of this inequality tends 
to zero, and (4.24) implies that the second term also tends to zero. 
Therefore, 

N{v{Ck),Ck)^N{v{Co),Co) inX*. (4.28) 

It follows from (4.10) that 

B*p{Ck)-B*p{Co) = N{v{Ck),Ck)-N{v{Co),Co) inX*. (4.29) 

Since the operator B* is an isomorphism from L2{^) onto V^, we obtain from (4.28) and 
(4.29) that 

p{Ck)^p{Co) in ^2(0), 

and the theorem is proved. ■ 
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We assign a mapping G : Hq{Q,) H ^{fl) as follows: 

H'om 9 C - G{0 = A{v{c),c) G H-\n). 

Here A : X x i?o(ri) —>■ H~^{^l) is given by (3.11) and v{() is the solution of the problem 
(4.9)-(4.11). 

Lemma 4.3. Suppose the conditions (CI), (C2), (2.7) -(2.9) are satisfied. Then G is a 
compact mapping of Hq{Q) into H~^{D,), the condition Ck Co in Hq{Q) implies G(Cfc) — >■ 
G(Co) in H-^{n). 

Proof. Let (^k — ^ Co in Hq{D,). Then Theorem 4.1 and the embedding theorem yield 

v{Ck)^v{Co) inUinr, (4.30) 
^ij (v iCk)) ^ Sijiv (Co)) inL2(r2), i,j = l,...,n, (4.31) 
Cfc^Co inUin). (4.32) 

In view of (4.31), (4.32) a subsequence {Crj} can be extracted from the sequence {Cjk} such 
that 

^ Co a.e. in (4.33) 
iHCr,)) ^ iHCo)) a.e. in n. (4.34) 

We represent the mapping G as a sum of two mappings G = Gi + G2, which we define as 
follows: 

(Gi(C),C)= [ {f + 0iui + vi0i)^dx-x7r{f,0, C^^H^i^), (4.35) 
(G2(C),C) = 2£ / [e{\Elf + C,x){X + I{u + v{0)-"^ 

+ V'l {I{u + v{0), \E\,f + C, x)] I{P{u + v{Q)) C dx, C, C e ^0 (^^) 

(4.36) 

If we will argue that 

Gi(Cik) ^ Gi(Co) in/f-i(0), (4.37) 

G2(Cik) ^ G2(Co) in/f-^(J^), (4.38) 

then the lemma will be proved. 
We denote 

hi = (r + Ck){ui + v{Ck)i) - (r + Co)(^ii + 'y(Co)i)- (4.39) 
It follows from (4.35) that 

(Gi(C,)-Gi(Co), i)\ = \l fki^,dx\ <{J2 [ fiid^yUlh' (4-40) 

1=1 

II • 111 being the norm defined by (3.9). 
(4.30) and (4.32) yield 

fki^O in L2(r2) as A; — 00, i = l,...,n. (4.41) 

By (4.40) and (4.41), we obtain (4.37). 
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By applying (3.1), (3.2), the Holder inequality, and the notations (4.18), we obtain 
^ (G2 (Cfe) - G2 (Co) , 1 = I _^ [Vkk I{P{u + v^)) 
-^00 1{P{u + v^'m dx = [ [^kk{P4 - Pe%){Pe% + ^4) 
+{^kk- ^m)I{P{u + v^))]idx < {Saik + a2k)U\\Le{n) < c((5Q;ijk + Q;2fe)||CI|i> 



Here 



(5 = ao A 2+02, 



aife = Yl (11^(4 -4)1^12(^)11^(4 +4)iUi2(n)> 



a2k = [ j {^kk - ^oo) 5 {IiP{u + v^))) 5 dx) ' 



(4.31) yields 



limaifc = 0. 



(4.42) 

(4.43) 
(4.44) 

(4.45) 
(4.46) 



By (4.33), (4.34) and the Lebesgue theorem, we obtain that limQ;2,j = 0. By analogy with 
the stated above, we can see that from any subsequence {Cm} extracted from the sequence 
{Cfe}, one can extract a subsequence {Ci} such that limQ:2i = 0. Therefore, 

lim a2k = 0. (4.47) 

(4.42), (4.46) and (4.47) imply (4.38), and the lemma is proved. ■ 

5. Proof of the Theorem 3.1. 
Consider the problem: find a function C such that 

c&H^m, (5.1) 

n{c,0--iM<C)X),0=o, ^eH^in), (5.2) 

X 

where v{Q is the solution of the problem (4.13). In this case there exists a unique function 
p{Q G L2(^) such that the pair {v{Q, p{C)) is the solution of the problem (4.9)-(4.11). 

It is obvious that, if C is a solution of the problem (5.1), (5.2), then the triple {v = 
^(0) P = P(C)) C) is a solution of the problem (3.15)-(3.18). 

Let |-2^A;| be a sequence of finite dimensional subspaces in HQiJl), such that 



2k CZ, 



fc+i) 



lim inf \\w — h\\i = 0, 

fc— >oo hGZf^ 



We seek an approximate solution of the problem (5.1), (5.2) of the form 

Ck^^k, 7r(Cik,0--(^(^(a),Cfc),6 = 0, ^eZk, 

where v{Ck) is the solution of the problem (4.13) for ^ = (^j^. 



(5.3) 



(5.4) 
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Applying the integration by parts, we obtain 

/ C{ui + Vi)—dx = - —{ui + Vi)Cdx = 0, 
Jn oxi Jq dxi 

veV, C^Hl{Q). (5.5) 
Taking into account (2.7), (2.8), (3.1), (3.2), (4.12) and (5.5), we reduce from (3.11) that 

-{a{v{Q,Q,(^) I < cIlClli, C G Hli^l). (5.6) 



Therefore, 



y(C) = - \{A{v{cuU) > IICII? - cllClli, (5.7) 

A. 



andy(C)>Ofor ||C||i > c. 

Prom the corollary of Brouwer's fixed point theorem [2] it follows that there exists a 
solution of the problem (5.4) and ||Cfc||i < c. Consequently, we can extract a subsequence 
{C^} from the sequence {Cfc} such that 

- Co in Hl{^). (5.8) 

Let r/o be a fixed positive number and ^ G ^^^y. By (5.8) and Lemma 4.3, we pass to the 
limit in (5.4) with k replaced by r]. Then, we get 

Co^HiiO), 7r(Co,e)--(^KCo),Co),6 = o, ee^^o. (5.9) 

(5.3) and (5.9) imply that the function ^ = Co is a solution of the problem (5.1), (5.2). The 
theorem is proved. ■ 

6. Generalized solution of the coupled problem for the function ip2. 
We denote 

1,- 

and determine a mapping Ai : X x Wq^{Q.) — ^ Q* as follows: 

{A^(y,C,U) = - [ ^^^-t^{ui + Vi)^dx 

Jn oxi 

+2e / [e{\E\,f + C,x){\ + I{u + v))-^2 
Jn 

+V'i(/(u + v), \E\,f + C, x)] I{u + v))i dx - X7r(f, 0, 

v^x, CeWo'^(O), CeQ. (6.1) 

1 — 6 

Here W^'^ (J^) is the closure of D{rt) for the norm of the Sobolev space 
The expression 



i=l 

defines a norm in VFq'^ (O) which is equivalent to the norm of W^^i{^). 
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Consider the problem: find a triple {v,pX) such that 

{v,p,OeX xL2{n)xW^'Hn), (6.3) 
{N{v,C),h)-{B*p,h) = {K + F,h), hex, (6.4) 
{Bv,q)=0, qeL2{n), (6.5) 

n{C,O--{Ai,{v,O,O=0, eeQ. (6.6) 
X 

The triple {u = u + v , p, t = f + () , where {v,p, ^) is a solution of the problem (6.3)-(6.6), 
is a generalized solution of the problem (2.1)-(2.6). 

Indeed, by use of the Green formula, one can verify that, if {v,p,Q is a solution of the 
problem (6.3)~(6.6), then u,p, r with u = u + v, and r = f + ^ is a solution of the problem 
(2.1)-(2.6) in the distribution sense. On the contrary, if u,p,T is a solution of the problem 
(2.1)-(2.6) such that (6.3) holds with v = u — u, (^ = t — f, then {v,p, Q is a solution of the 
problem (6.3)-(6.6). 

Theorem 6.1. Let be a bounded domain in R", n = 2 or 3 with a boundary of the class 
C^. Suppose the conditions (CI), (C2), (2.7)- (2.9) are satisfied. 
Assume also that 



inf — ||n + w^llifira-)" + f ll-fC + Flly* + 2(aoA 2 +a2)f / Iiu + w)dx\^\ 
wev X L ^ ' 2ai\ \Jq J J 



<1, 
(6.7) 

where c is a constant depending on the domain Q. Then, there exists a solution of the 
problem (6.3) -(6.6) 

Proof. 

l)Lct {Pk}kLi be a sequence of regularizing operators assigned by 



where 



Pku{x) = / u>k{\x — x'\)u{x')dx' , X € ri, (6-8) 

Wjfc G C°°(R+),suppa;fe = [0,6jk], ujk{z)>0, 2; G M+, / Uk{\x\)dx = 1, 

•^K" (6.9) 

and also lim bj- = 0. In this case we have 

lim\\PkU — u\\fji(^Q-jn = 0, uG-ff^(r2)", 

lim\\PkU-u\\L^(^a)n = 0, ueL2{nT, (6.10) 

and it is assumed here that the function u is prolonged in J7i D ri, such that u belongs 
to H^{Q,iY or L2{^iY, respectively. 

By virtue of the Theorem 3.1 for each operator P = Pj. there exists a solution of the 
problem (3.15)-(3.18), which we denote by v'', pk, Cfe- In this case the function is a 
solution of the following problem: 

aeH'om, Y.^ = f^'^-Y.Q^2 ^^D*{n), (e.n) 

1=1 * i=l * 



where 



Pk = -i^i + vf f^''^'^ - - Vkk miu + v'')), (6.12) 
X 9xi X 
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and (fkh is defined in (4.18). 
We have 



[Ui + Vi 



where 



M 



Ik 



df_ 
dxi 



dxi 



Li(n) 



Li(n) 



< Mife + M2k, 



2k 



(6.13) 
(6.14) 



It follows from (2.7)-(2.9), (4.4) and (4.12) that 

df 



i=l 



Mik = ^\\ui + v^\\Le(n) < ci||?i + / 11^1(^2)" ll^lb < C2- (6.15) 



(2.7), (4.4) and (4.12) imply 



dCk 



1=1 



M2k < XI \\ui + Vi\\Leia) r.(a\- ^adl^^lkM^^)" + lb^lU)IICfe||2 < CscllCfclb,^ 



L6(n) 

5 



where 



2a^ 



\K + F\\v* +2(aoA-2 +a2)( j I 



{u)dx^ . 



Taking note of (6.13), (6.15) and (6.16), we obtain 

.M d{f + Ck 



{Ui + <) 



Li(n) < C3c||Cfe||2 + C2. 



(6.16) 
(6.17) 

(6.18) 
(6.19) 
(6.20) 



(2.7), (3.1), (3.2) and (4.12) yield 

llm I{Pk{u + ^^''))||Li(n) < C4. 
It follows from (6.12), (6.18) and (6.19) that 

||A||Li(n)< ^Iiail2 + C5. 
X 

Since the boundary S of the class and 

Q2f J. 
-^^H'^{a)^W-^^i{Q.), i = l,...,n, 

where VF^^'t(f7) is the space dual to Wq^'^(J7), we obtain from (6.11) and known results 
[13], [14], Section 1.1 that 

||ai|2<C6(||A||^_,a^^^ + ||f||^,,a^^^). (6.21) 

The embedding of Li(fi) into TV'^^'^(ll) is continuous because Wq^'^(O) C C{Vt) for n = 2 
and 3. So that by (6.20) and (6.21), we obtain 



Kkh < c7\\l3k\\Li{n) + C8 < — llCfclb + cio- 



That is 



(i-^)iiaii2<cio 



(6.22) 



(6.23) 



It is evident that, if a function u satisfies the conditions (2.7), then the function u + w, 
where w G V, satisfies (2.7) also. Therefore, we can consider that, in the operators A^, A and 
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Ai the function u is replaced by the function u + w, where u; is a function from V. In this 
case, if {v,p, C) is a solution of the problem (6.3)-(6.6) for the function u, then (v — w, p, Q 
is a solution of (6.3)-(6.6) for the function u + w, and instead of (6.23), we obtain the 
following inequality 

[l-7H]||ai|2 <cio, (6.24) 

where, see (6.17), 

c 1 1 f ~ 

j{w) = — \\u + w\\Hi(Q)n + - — (\\K + F\\v* +2{aoX''^ + a2)( / I{u + w)dx)^) . 

X L ' ' 2ai\ \Jn ^ ^(^ 25) 

If the condition (6.7) with c = cg is satisfied, then there exists a function w e V such that 

7(«;) < 1, and by virtue of (6.24) the sequence {("fc} is bounded in Wq'''{Q). 

We consider that the function u is replaced by the function u + w in the operators N, A 
and Ai, and we still denote by u the function u + w. Then for the new function u, we have 
in (6.23) cgcx~^ < 1- Therefore, a subsequence {v"^, Pm, Cm} can be extracted from the 
sequence {v'', Pk, Cfc} such that 

Cn,^Co inW^'^n), (6.26) 
Cm Co in a.e. in Q, (6.27) 

ym _^ ^0 -^^ 28) 

N{v"',Cm)^cx inX*, (6.29) 

Pm^Po inL2{ft), (6.30) 

2). Now we are concerned with the passage to the limit. It follows from (6.4) and (6.5) 
that 

{N{v"',Cm),h)-{B*pm,h) = {K + F,h), hex, (6.31) 
{Bv"',q)=0, q£L2in). (6.32) 

Observing (6.28)-(6.30), we pass to the limit in (6.31), (6.32), and obtain 

a-B*po = K + F inX*, (6.33) 
div'y° = 0. (6.34) 

Lemma 4.1 implies 

{N{v"',Cm)-N{g,Cm),v"'-g)>0, g e X. (6.35) 
Taking into account that (pjn,Bvjn) = 0, by (6.28) and (6.31), we obtain 

{N{v^, Cm),vn = {K + F, v^) ^{K + F, v""), (6.36) 

and 

\im{N{v^,Cm),9)-{B*po,g) = {K + F,g), g e X. (6.37) 

By using (CI), (C2), (6.27) and the Lebesgue theorem, we deduce 

[e(|£;|, f + Cm, x){\ + I{u + 5))-5 + V'i(/(n + g),\E\,f + Cm, x)] 

xeij{u + g) [e{\E\,f + Co, x){\ + I{u + g))-^ 
+il^i{I{u + g),\E\,f + Co,x)]eij{u + g) in L2(J^) as m ^ oo. (6.38) 
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(6.28) yields 

mL2(0). (6.39) 

It follows from (3.10), (6.38) and (6.39) that 

(7V(ff,C™),^™ - <?) ^ (iV(5,Co),^° - g). (6.40) 

Observing (6.36), (6.37) and (6.40), we pass to the limit in (6.35), this gives 

(^K + F-N{gXo)+B*p^,v'' -g)>Q, g e X. (6.41) 

We choose g = v^ — 'yh, hGX,j>0, and consider 7^0. Then, (4.1) implies 

iK + F-N{v^,Co) + B*po,h)>0. (6.42) 

This inequality holds for any h & X. Replacing hhy —h shows that here the equality holds 
true. Consequently, the triple {v = , p = po, ( = Co) meets the equations (6.4) and (6.5). 
It follows from (4.2) that 

{Niv^, Cm) - N{v'>,Cm),v^ - v"") > fi2\\v'^ ' v'^fx- (6.43) 

Granting (6.36), (6.37), (6.40) one can recognize that the left-hand side of (6.43) tends to 
zero. Therefore, 

and we can regard that 

liv"^) I{v^) a.e. in n. (6.45) 

Otherwise, we can extract a subsequence, still denoted by {v"^}, such that (6.45) holds true. 
3). The pair {v"^,(m) meets the following equation (see (3.11) and (3.18)) 

7r(Cm,0--(t^m,0 = 0, CeH^{n)nL^{n) = Q, (6.46) 

X 

where 

Jn (^^i Jn (-g 47") 

iprnm being defined in (4.18). 

By virtue of (6.44) v"" in L6(n)", so that by (6.26), we obtain 



lim / ^il±M (^i, + dx= [ ^^^±^ (u, + v^)^ dx, C G Q. 



We have 



where 



(6.48) 



(6.49) 



being defined in (4.18). 



MT= [ {fmm - V00)I{U + v'')dx , (6.50) 

Jn 

M'T = \ E V^mmiPmS^ ' e%){PmeTj + ^l)dx , (6.51) 
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Taking into consideration (CI), (C2), (6.27) and (6.45), we deduce from the Lebesgue 
theorem that 

UmX^ = 0. (6.52) 

(CI) and (C2) imply 

n 

< (aoA-5 + as) ^ - e^h^in) 

n 

X \\Pm£ij + £ij\\L2in) < C ^ WPm^ij - £ij\\L2{n) 

n 

< c J2 - 4)IU2(n) + WPmsl - 6° ||L,(n)) (6.53) 

It follows from (6.10) and the Banach-Steinhaus theorem, that the norms of the operators 
Pm are uniformly bounded in i?i(J))" and L2{n). So that (6.10), (6.44) and (6.53) yield 
limTW™ = 0. Therefore, the pair {v = v^, ( = (q) satisfies the equation (6.6), and the 
theorem is proved. ■ 

Theorem 6.1 is a generalization of the result obtained in [17] for a model of non-linear 
viscous fluid. 



7. Coupled problem for the viscosity function (pi. 

For the viscosity function tpi defined by (1.11), we take the equation of the balance of 
thermal energy in the following form: 

X^Q^ + '^eM\ElK^^,E),IiP^),r)I{Pu)-Ui— =0 in Q, (7.1) 

i=l * * 

where 

M\E\,^iiu, E),I{Pu),t) = b{\E\,^i{u, E),T)i\ + /(Pn))-^ + ^P{I{Pu), \Elfi{u, E),t), 

(7.2) 

and P is an operator of regularization determined in (3.12), (3.13). 

The presence of the operator P in (7.1) means that the model is not local and it is natural 
from the physical point of view (see Section 3). 

The motion equations, the condition of incompressibility and boundary conditions are 
given by (2.1) under k = 1, (2.2), and (2.4)-(2.6). 

Wc assume that (2. 7) -(2. 9) hold and the functions b and i/j satisfy the following conditions 
which are similar to the conditions (CI) and (C2). 

(Cla): b : {yi,y2, Vs) b{yi,y2, y^) is a function continuous in M_|_ x [0, 1] x M, and in 
addition 

< %i, y2, ys) < ao, (yi, 2/2, 2/3) e x [0, l] x R. (7.3) 

(C2a): V : (2/1,2/2,2/3,2/4) ^ V' (2/1, 2/2, 2/3, 2/4) is a function continuous in M+ x R+ x 
[0,1] X M, and for an arbitrary fixed (2/2,2/3,2/4) ^ 1^+ [0,1] 1^ the function 
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V'(-) 2/2; 2/4) : 2/1 ^^(2/1) 2/2) 2/3) 2/4) is continuously differentiable in M_|_, and the 
following inequalities hold 

a2 > ^{yi,y2, y3,yi) > oi, (7.4) 
dtp 

V'(2/i, 2/2,^3,2/4) + 2— (2/1,^2,2/3,^4)2/1 >a3, (7.5) 
oyi 

4^(2/1, 2/2, 2/3, 2/4) 2/1 < «4 (7.6) 
02/1 

We introduce mappings A^i : X x H^{n) X* and ^1 : X x H^{n) H''^{n) as follows: 

(iVi(^;,C),^) = 2 / [b{\E\,^i{u + v,E),f + 0{X + I{u + v)y'^ 
Jn 

+i/j{I{u + v), \E\,fx{u + v,E),f + C)]£ijiu + v)eij{h) dx, 

v,hex, CeH^in), (7.7) 

{A,{v, C), = - _^(«i + ^i)^^^^ ? 

+2e [ M\E\,fi{u + v,E),I{P{u + v)),f + OHPiu + v))^dx-XTr{f,0, 
Jn 

vex, c,eGi^o(^)- (7.8) 

Consider the following problem: find a triple {v,p,Q such that 

{v,p,(:)eXxL2{n)xH'o{n), (7.9) 

iNi{v,0,h)-iB*p,h) = iK + F,h), hex, (7.10) 

{Bv,q)=0, qeL2{p), (7.11) 

7r(C,0-^(^i(^,C),e) =0, ieHl{^). (7.12) 

Here we use the notations (3.19). 

If {v,p, C) is a solution of the problem (7.9)-(7.12), then {u + v,p,f + Q is a generalized 
solution of the problem (2.1) for k = l, (2.2), (2.4)-(2.6) and (7.1). 

Theorem 7.1. Let Q be a bounded domain in R", n = 2 or 3, with a boundary S of the 
class , and suppose the conditions (Cla), (C2a) and (2.7), (2.9) are satisfied. 
Assume also that there exists a function f such that (2.8) holds and in addition 

f e L^in). (7.13) 

Then exists a solution of the problem (7.9)-(7.12) 



8. Lemmas and estimation of the term generated by the convection. 

Below we suppose that is a bounded domain in M", n = 2 or 3 with a boundary S of 
the class C^. 
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Lemma 8.1. Let p{x) he the distance from x to S. Then for an arbitrary 7 > there exists 
a function G (O) such that 

a^y = 1 in some vicinity of S depending on 7, (8-1) 



a^{x) = at p{x) > 25(7), (^(7) = exp( ), 



d 



dxk 



< 



7 



p{x) 



7 



at p{x) < 25(7), k = 1, . . . ,n. 



(8.2) 
(8.3) 



The function satisfying the above conditions was constructed by Hopf [4], see also [7]. 
Lemma 8.2. There exists a constant c depending only on Q such that 

<c||ii;||i, weH^{n). (8.4) 



1 

- w 
P 



L2(,n) 



This lemma is proved by a partition of unity and local maps followed by the application 
of the Hardy inequality, see e.g. [18], Lemma 1.10, Chapter 2. 

Lemma 8.3. Suppose a function f satisfies (2.8) and (7.13). Then for an arbitrary /? > 
one can construct a function G such that 

GeH\n)nL^{n), G|^ = f|^, 

and in addition 



C— 
dxi 



</?IICI|i, C^Hl{^), ^ = l,...,n. 



(8.5) 
(8.6) 



Proof. Consider the function 



G = a^f. 



(8.7) 



This function satisfies (8.5). It follows from (7.13), (8.7) and Lemma 8.1 that 

df 



,n. 



(8.8) 



(8.8) implies 



^dG 
' dxi 



< ci 7 



< 02(^7 



Lein) 

"5 



+ 



+ 



P<2<5(7) 



. df 
' dxi 



' dx 



where 



IICI|L3(n)<^i(7)), 



df \2 



(8.9) 



^.10) 



i = 1, . . . , n. 

/p<25(7) ' ' 

Here we have applied the Holder inequality with the indexes | and | to the integral on the 
right-hand side of (8.9). 
(8.9) and Lemma 8.2 yield 



dG 
dxi 



<C3||CI|i(7 + f^i(7)), 



i = 1, . . . ,n. 



(8.11) 
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Since ^ G -Z^2(^), we obtain from (8.10) that <Tj(7) ^ as 7 ^ 0, and the Lemma is 
proved. ■ 

We consider a mapping A2 :V x Hq{Q,) H~^{Cl) that is defined as follows: 



{A2iv,0,0 = - |^{i^i + v^)^^^^^dx, vex, C,^&H^{^)- 



(8.12) 



Lemma 8.4. Suppose the conditions (2.7), (2.8) and (7.13) are satisfied. Then for an 
arbitrary rj > one can determine the function G so that (8.5) is satisfied, and in addition 

|(^2(t',C),C)l<c||C||i+r?(||^;||^ + ||C||?), veV, C^H^i^)- (8-13) 

Proof. It follows from (8.12) that 

{A2{v,0,0 = Mi + M2 + Ms, (8.14) 

where 



Ml 



Ui—Cdx, 



M2 



(2.7) and (8.6) yield 



\Mi \ <^\\ui\\L^(u) 



i=l 
n 



\M2\<J2MLe{^) 



i=l 



dG ^ 

dxi 

dG ^ 
dxi 



Lain) 



< /3ci||u||ii-i(n)n||CI|i < C2||CI|i> 



Since div('u + v) = 0, we obtain by integration by parts that 



Ms = - / ( 

Jo, 

Now (8.13) follows from (8.16)-(8.18) and Lemma 8.3. 



OXi Jn 



dxi 



dx = 0. 



(8.15) 

(8.16) 

(8.17) 
(8.18) 



9. Proof of Theorem 7.1. 

1) It follows from the proof of Theorem 4.1, that if {v,p,Q is a solution of the problem 
(7.9)-(7.12), then 



l^^llx < —{\\K + F\\v'+ii4). 



(9.1) 



By (3.12), (3.13) and (9.1), we conclude that there exists a constant 61 such that 

mP{^ + vmcia)<bi- (9.2) 

Consider a function 993 such that (^3 is a function continuous and bounded in R+ x [0, 1] x 
Mi X M — s- Ml and in addition 



(pi{z) at z = {zi,Z2, Z3, za) G M+ X [0, 1] x [0, 61] x : 
at 2; = (zi, 2:2, 2^3, Zi) G M+ X [0, 1] X (62, 00) X M, 



(9.3) 
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where 62 > 5i and ipi is the function defined by (7.2), i.e. 

ipi{z) = b{zi,Z2, Z4){X + Z3)"3 + ^(Z3, Zi, Z2, Z4). (9.4) 

Define operators N2 : X x H^{n) X* and A3 : X x H^{n) H-^{n) by 
{N2(,v,0,h) = 2 [ [b{\E\,fi{u + v,E),G + 0{X + I{u + v)r'^ 

+'iIj{I{u + v), \E\,iJ,{u + v,E),G + C)]£ij{u + v)eij{h) dx, 

v,heX, C&H^i^), (9.5) 

(^3(^;, C), = 2£ / M\E\,K^ + V, E), I{P{u + v)), G + Ol{P{u + v))i dx - x7r(G, 0, 

i;ex, c,Ce^^o(^^), (9-6) 

where G is the function defined in Lemma 8.3. 

We consider the problem: find a triple of functions {v,p,9) satisfying 

{v,p,e) e X X L2{n) X H^{n), (9.7) 

{N2{v,e),h)-{B*p,h) = {K + F,h), hex, (9.8) 
{Bv,q)=0, q£L2m, (9.9) 

7T{e,o - -[iMv,o),0 + {Mv,&),0] =0, C e H^{n). (9.10) 

X 

It follows from the proof of theorem 4.1 that the inequality (9.1) is also fulfilled for the 

solution of the problem (9.7)-(9.10). Therefore, (9.2) holds, and if {v,p,9) is a solution of 
the problem (9.7)-(9.10), then the triple {v,p,C, = G — f + 9) is & solution of the problem 
(7.9)-(7.12). 

The problem (9.7)-(9.10) is equivalent to the following problem: find a pair {v,9) such 
that 

{v,9) eV X Hl{Vl), (9.11) 
{N2{v,9),h) = {K + F,h), heV, (9.12) 

7r{9,0 - -[{Mv,d),0 + {Mv,d),0] =0, C e H^i^). (9.13) 
X 

Indeed, if {v,p, 9) is a solution of the problem (9.7)-(9.10), then the pair (d, 9) is a solution of 
the problem (9.11) (9.13). Conversely, let {v,9) be a solution of the problem (9.11) (9.13). 
By virtue of (9.12) and Lemma 4.2 there exists a unique function p G 1/2 (i^) such that the 
triple {v,p,9) is a solution of the problem (9.7)-(9.10). 

2) Let us prove the existence of a solution of the problem (9.11)-(9.13). 

We define a mapping U:Vx H^{n) ^ V* x H-'^{n) by 

{U{v, 9), {h, 0) = iN2{v, 9),h)+ 7r{9, - -[{M^, 0),O + (Mv, 0),O] - {K + F, h), 

X 

v,heV, 9,CeH^{n). (9.14) 

It is easy to verify that the problem (9.11)-(9.13) is equivalent to the following one: find a 
pair (v, 9) such that 

{v,9) eV X Hl{Vl), (9.15) 
{U{v,9),{Ki))=Q, {h,OeVxHl{n). (9.16) 
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Indeed, taking the pair (/i, 0) in (9.16), we obtain (9.12), and taking the pair (0,^), we get 
(9.13). Conversely, by adding (9.12) and (9.13), we get (9.16). 

Let {Vm} and {Z^} be sequences of finite-dimensional subspaces in V and Hq{^1), re- 
spectively, such that 

lim inf \\u - g\\x = 0, ueV, (9-17) 

m— »oo gGVm 

lim inf = 0, w e H}(n), (9.18) 

Vm C Vm+1, Zra G ^m+L (9.19) 

We search for an approximate solution v^, of the problem (9.11)-(9.13) in the form 

{Vm. dm) eVmX Zm, (9.20) 

{N2ivm,9m),h) = iK + F,h), heVm, (9.21) 

),0 + (^3( ),0]=0, i^Zm. (9.22) 

X 

(9.3) and (9.6) imply 

|(^3(^^,^),^)| <ci / \e\dx + x{AG,G))He\u<c2\\e\\^. (9.23) 

Inequality (4.3) is also realized for the operator N2. So that by using (8.13), (9.14) and 
(9.23), we obtain 

{u{v,e), {v, e)) > + \\o\\l - calif lU - C4||0||i - -(11^^113^ + II^IlD, 

X 

{v,e) &V Hl{n). (9.24) 

By Lemma 8.4 we can reckon that 

1 

V = 2Xmm(/X3,l). 
Then there exists a constants r > such that 

{U{v,e),{v,e))>Q, if ||i;|U + ||^||i >r. (9.25) 
Therefore, for an arbitrary m G N there exists a pair (?;^, 0^) satisfying 

),(/i,e))=0, {h,i)eVm>iZm. (9.26) 
This pair is the solution of the problem (9.20)-(9.22), and moreover 

Ibmllx + ll^mlli <r. (9.27) 
It follows from the proof of Theorem 4.1 that 

\\vm\\x<—{\\K + F\\v*+m). (9.28) 
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3) By virtue of (9.27), we can extract a subsequence {vh, 0^} from the sequence {vm, 9m} 
such that 

Vk^vo in V, (9.29) 

Vk — ^0 in L4{n)'" and a.e. in CI, (9.30) 

9k in H^{n), (9.31) 

9k 9o in L2{^) and a.e. in Q, (9.32) 

iV2(^;ifc,^ifc)^JP- in X*. (9.33) 

Let /cq be a fixed positive number, and h € Vk^. Observing (9.33), we pass to the Umit in 
(9.21) with m replaced by k, and obtain 

{J^,h) = {K + F,h), heVk,. (9.34) 

Since ko is an arbitrary positive integer, we get by (9.17) and (9.34) that 

J^ = K + F in V*. (9.35) 
We present the operator N2 in the form 

N2(,v,0 = N{v,v,0, (9.36) 

where the operator {v, w, Q — N{v, w, Q is considered as a mapping ofV xV x Hq{Q,) into 
V* according to 

{N{v, w, 0, h) = 2 f [b(\E\,fi{n + w, E), G + C)(A + I(u + v))-'^ 

+ip{I{u + v), \E\,ii{u + w, E),G + C,)]£ij{u + v)eij{h) dx, 

v,w€V, CeH^{n). (9.37) 

For an arbitrary fixed w & V and ( G Hq{Q,) the operator N{.,w,z) : v N(v,w,z) is 
monotone (see (4.2)). Therefore, 

{N{vk,Vk,ek)-Niz,Vk,ek),Vk-z)>0, zeV, fcGN. (9.38) 

(9.30), (9.32) and Lebesgue theorem give 

lim\\N{z,Vk,9k)-N{z,vo,9o\\x* =0, z e V. (9.39) 

By (9.21), (9.29) and (9.36), we obtain 

{Nivk, Vk, 9k), Vk) = {K + F, Vk) ^{K + F, vo), (9.40) 

and (9.17) yields 

hm{N{vk, Vk, 9k), z) = {K + F, z), zeV. (9.41) 

Observing (9.29), (9.39)-(9.41), we pass to the hmit in (9.38). 
Then we get 

{K + F-N{z,vo,9o),vo-z)>0, z € V. (9.42) 

We choose z = vq — 7/1, 7 > 0, /i G and consider 7^0. Then, taking into account that 
A^(., Vq, 9o) : z — > N{z, vo, 9q) is a Lipschitz continuous mapping of V into V* (see (4.1)), we 
obtain 

{K + F-N{vo,VQ,9Q),h)>Q. (9.43) 
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This inequality holds for any h eV. Therefore, replacing h by —h shows that equality holds 
true in (9.43). Prom here by (9.29), (9.31) and (9.36), we deduce 

{vo,eo)eVxH^{n), {N2{vo,eo),h) = {K + F,h), heV. (9.44) 

(8.12), (9.6), (9.30)-(9.32) imply 

lim(^2(^^ifc,^ifc),6 = {A2{vo,eo),0, 

lim(^3(«ik, dk), = (Mvo, do), 0, C e H^m. (9.45) 

Bearing in mind (9.31) and (9.45), we pass to the limit in (9.22) for fixed ^ G Zk^^. Then by 
(9.18), we deduce 

A0o,O - -[{Mvo,do),0 + (^3(^0,^0), 6] = 0, C G H^{n). (9.46) 

It follows from (9.44) and (9.46) that the pair (v = vq, 9 = 6q) is a solution of the problem 
(9.11)-(9.13). Hence, there exists a unique function p such that the triple {v = vq, p, 9 = 9o) 
is a solution of the problem (9.7)-(9.10). 
(9.28) and (9.29) yield 

HWx <—{\\K + F\\v*+fi4). (9.47) 

Since (9.1) implies (9.2), we get \\I{P{u + vo))\\c(Ti) < bi- From here by (7.8), (9.3) and 
(9.6), we obtain that the triple (;; = vq, p, = G — f + 9o) is a solution of the problem 
(7.9)-(7.12). The theorem is proved. ■ 
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